Spin polarization decay in spin- 1/2 and spin-3/2 systems 
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We present a general unifying theory for spin polarization decay due to the interplay of spin 
precession and momentum scattering that is applicable to both spin-1/2 electrons and spin-3/2 holes. 
Our theory allows us to identify and characterize a wide range of qualitatively different regimes. For 
strong momentum scattering or slow spin precession we recover the D'yakonov-Perel result, according 
to which the spin relaxation time is inversely proportional to the momentum relaxation time. On the 
other hand, we find that, in the ballistic regime the carrier spin polarization shows a very different 
qualitative behavior. In systems with isotropic spin splitting the spin polarization can oscillate 
indefinitely, while in systems with anisotropic spin splitting the spin polarization is reduced by spin 
dephasing, which is non-exponential and may result in an incomplete decay of the spin polarization. 
For weak momentum scattering or fast spin precession, the oscillations or non-exponential spin 
dephasing are modulated by an exponential envelope proportional to the momentum relaxation time. 
Nevertheless, even in this case in certain systems a fraction of the spin polarization may survive 
at long times. Finally it is shown that, despite the qualitatively different nature of spin precession 
in the valence band, spin polarization decay in spin-3/2 hole systems has many similarities to its 
counterpart in spin-1/2 electron systems. 
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I. INTRODUCTION 

The achievement of a lasting spin polarization has 
been a long-standing goal in semiconductor physics. Suc- 
cessful efforts to generate a spin polarization magnet- 
ically, optically and electrically, have yielded a steady 
stream of novel physics and promising applications. 1,2 
Ferromagnetic semiconductors are edging towards room 
temperature^ and spin currents have been measured 
directly. 4 Successes such as these have turned semicon- 
ductor spin electronics into a vibrant and rewarding area 
of research, as well as a promising candidate for novel 
information processing methods. 

Both for fundamental physics and for technological ap- 
plications, it is important to know how to maintain a spin 
polarization once it is generated. Therefore, a detailed 
understanding of the mechanisms leading to spin polar- 
ization decay is critical in all areas mentioned above. In 
the return to equilibrium of an excess spin polarization 
spin-orbit interactions play an important role. Spin-orbit 
coupling always gives rise to spin precession, and the in- 
terplay of spin precession and momentum scattering is 
frequently the main cause of spin polarization decayi^* 6 -! 7 - 
A spin polarization in a semiconductor may also decay 
via spin flips induced by momentum scattering or by 
exchange interactions, though these mechanisms have a 
more limited range of applicability ^ 8 ' 9 ' 10 

Spin relaxation in spin-1/2 elec- 
tron systems has received considerable 

at |- CI1 ^ on 5.6.7.8.9 .10.11.12.13.14.15.16.17.18.19.20.21.22.23.24.25.26 .^ 

For electrons the spin-orbit interaction can al- 
ways be represented by a Zeeman-like Hamiltonian 
H = (fi/2) a ■ O(fe) describing the interaction of the 
spin with an effective wave vector-dependent magnetic 
field fi(fc). The electron spin precesses about this field 
with frequency Q = |0(fe)|. An important parameter is 



the product of the frequency Q times the momentum 
relaxation time t p . In the ballistic (clean) regime 
no scattering occurs and the temperature tends to 
absolute zero, so that f2r p — ■> oo. The weak scattering 
regime is characterized by fast spin precession and little 
momentum scattering due to, e.g., a slight increase 
in temperature, yielding flr p 3> 1. In the strong 
momentum scattering regime flr p <C 1. 

Electron systems are often in the strong scattering 
regime. In this case the main mechanism leading to 
spin polarization decay is the D'yakonov-Perel (DP) 
mechanism ^£2. which was shown to be dominant over a 
wide range of temperatures 7 - and, for particular forms of 
fi(fc), to lead to a noticeable anisotropy in the relaxation 
times for different spin component s 13 ' 14 and anisotropic 
spin diffusion*^ Most past work has concentrated on 
this regime. On the other hand, in recent years state- 
of-the-art technology has enabled the growth of ballis- 
tic samples which have been at the forefront of spin- 
related experiments* 1 ^ Yet spin polarization decay in 
ballistic spin-1/2 systems has received comparably little 
attention 1 - ' 20 ' 21 and has been considered recently mostly 
in the context of spin transport in an electric field i 17 ' 18 

For spin-3/2 holes the spin-orbit interaction cannot be 
written as an effective field, and spin precession is quali- 
tatively different j^ 9 - Since spin-orbit coupling is more im- 
portant in the valence band, hole spin information is lost 
faster, and the relative strengths of spin-orbit coupling 
and momentum scattering can vary. Yet spin relaxation 
2©f spin-3/2 holes has also been studied to a lesser extent, 
both experimentally^ - and theoreticallyi& 3 i 3 ^i 3 ^ 4 .i2 5 . A 
theory of spin relaxation valid for electrons and holes 
in all regimes of momentum scattering does not, to our 
knowledge, exist to date. 

With these observations in mind, we present in this 
article a general unifying quantitative theory for the re- 
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turn to equilibrium of excess spin polarizations in the 
conduction and valence bands of semiconductors brought 
about by the interplay of spin precession and momentum 
scattering. We do not rely on the assumption, made in 
most previous work,M i Z & . :u i 12 i 13 1 14 i 15 i 16 i 17 l that !1t p «1. 
We demonstrate that spin polarization decay in differ- 
ent regimes of momentum scattering in spin-1/2 electron 
and spin-3/2 hole systems contains considerable rich and 
novel physics. For example, spin polarization decay has 
often been assumed to be proportional to e~ l l Ta , where r a 
is referred to as the spin relaxation time. However, if the 
magnitude of the spin-orbit interaction is anisotropic (as 
is usually the case in systems studied experimentally), 
spin-polarization decay can occur even in the absence 
of momentum scattering. This process is characterized 
by a non-exponential decay and is sensitive to the ini- 
tial conditions, and cannot therefore be described by a 
spin relaxation time. Weak momentum scattering intro- 
duces a spin relaxation time t s cx r p (unlike strong mo- 
mentum scattering, which gives the well-knowr&l trend 
t s oc t~ ) , yet even in the presence of weak momen- 
tum scattering a fraction of the polarization may sur- 
vive at long times. It will emerge from our work that, 
in the ballistic and weak momentum scattering regimes, 
the concept of a spin relaxation time is of very limited 
applicability and in general does not provide an accurate 
description of the physics of spin polarization decay. 

We emphasize that the results presented in this paper 
are true for (delocalized) electron spins in any nonmag- 
netic solid where spin-orbit coupling is important. Since 
in today's experiments mobilities range over many orders 
of magnitude, the results presented arc directly relevant 
to ongoing state-of-the-art research. 

The outline of this article is as follows. In section II we 
discuss the time evolution of the density matrix, deriving 
an equation which describes the return to equilibrium of 
a spin polarization. We demonstrate that in the general 
case there exists a fraction of the spin polarization which 
does not precess, and explain its relevance to the sub- 
sequent time evolution of the spin polarization. Section 
III is devoted to spin-1/2 electron systems, in which first 
the known D'yakonov-Perel' limit is discussed, then the 
complex situations in the ballistic and weak momentum 
scattering regimes are presented. We stress the impor- 
tance of non-exponential decay and of incomplete spin 
dephasing. Finally, in the last part we demonstrate that, 
although spin precession is qualitatively different in spin- 
1/2 electron and spin-3/2 hole systems, spin polarization 
decay in these systems can be understood based on the 
same fundamental concepts. 



II. TIME EVOLUTION OF THE DENSITY 
MATRIX 

We assume a nonequilibrium spin polarization has 
been generated in a homogeneous, unstructured system 
and study its time evolution in the absence of exter- 



nal fields. The system is described by a density matrix, 
which in principle has matrix elements diagonal and off- 
diagonal in momentum space. Since the spin operator 
is diagonal in the wave vector fc, we will only be con- 
cerned with the part of the density matrix diagonal in 
momentum space, which is denoted by p. Henceforth, by 
"density matrix" we understand the part of the density 
matrix diagonal in wave vector. 

The spin density is given by (S) = trSp — trSp, 
where S is the spin operator, and the overlinc represents 
averaging over directions in momentum space. Only the 
isotropic part p of the density matrix is responsible for 
spin population decay^ It is therefore convenient to di- 
vide p into p = p + g, where g is the anisotropic part of 
p. Based on the quantum Liouville equation, we obtain 
an equation describing the time evolution of p (Ref. [8|), 
which in turn is split into a set of equations for p and g 
similar to those found by Pikus and Titkov>£ 
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These equations hold both for spin-1/2 electrons and for 
spin-3/2 holes. We assume elastic scattering by short- 
range impurities, implying that the collision term involv- 
ing p vanishes^ and the remainder is proportional to the 
inverse of the scalar momentum relaxation time^> 1/t p . 

Before proceeding, we would like to make two remarks 
concerning the form of the scattering term. Firstly, in 
the presence of spin-orbit coupling both intraband and 
interband transitions exist, while we have assumed a sim- 
plified form of the scattering term. In the version of 
the relaxation time approximation employed in this work 
the spin splitting of the bands is not taken into account 
in the scattering term. This approximation is justified 
by the fact that spin eigenstates are generally not en- 
ergy eigenstates, and it can be straightforwardly shown, 
based on the theory we present, that accounting explic- 
itly for interband transitions will not change the funda- 
mental physics of spin polarization decay, rather it will 
only give less transparent solutions. Furthermore, spin- 
flip scattering in nonmagnetic systems is third-order in 
the scattering potential and/or first order in the ratio of 
the spin-orbit splitting and the kinetic energy 

Secondly, it should be noted that, for degenerate carri- 
ers, the return to equilibrium requires energy dissipation. 
However, as noted above, in a nonmagnetic material 
with spin-orbit coupling the spin eigenstates character- 
izing the nonzero spin polarization are not energy eigen- 
states. On the other hand, unlike in, e.g., nuclear sys- 
tems, the nonthermal energy characterizing this nonequi- 
librium configuration is essentially a kinetic energy, but 
it is not in the spin degree of freedom. Therefore, energy 
dissipation has no qualitative effect for the main conclu- 
sions in our paper. 

A solution to Eq. (|lb[) can be obtained by making the 
transformation g = e~ lHt / h gi e lHt / h , which is analogous 
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to the customary switch to the interaction picture. This 
transformation turns Eq. (|lb[) into an equation for gj 

at t p ~ at' 1 > 

where pi is defined by p = e~ %m l h pi e lHt / h . Treating the 
RHS as a source term, this equation allows an analytical 
solution using an integrating factor. Substituting this 
solution into Eq. ([Taj) yields 

+ [ dtfe-Q-*'^' e-^*-*')^ [H, p{t')] e^(*-*')A 
at hT p J 

= -1 e -*/r P e -iHt/h \ H po l e iHt/h 

h 

(3) 

where po is the initial value p(t = 0). This equation is 
the main result of our paper. It describes the precession- 
induced decay of spin polarization in all regimes of mo- 
mentum scattering for any nonmagnetic solid state sys- 
tem with spin-orbit interactions. This equation does not 
anticipate any particular form of spin polarization decay, 
such as exponential decay. 

The form of the initial density matrix po is impor- 
tant and lies at the root of the novel physics discussed 
in this paper. In general po has two contributions, 
Po = Po\\ + Po_l- The component poll commutes with H 
and is given by p || = ( tr poH/tr H 2 ) H, in a generaliza- 
tion of Gram-Schmidt orthogonalization. poj_ is simply 
the remainder, and it satisfies the condition tr po±H = 0. 
Po|| is a matrix that is parallel to the Hamiltonian, and 
represents the fraction of the initial spin polarization that 
does not precess, or alternatively the fraction of the ini- 
tial spins that are in eigenstates of the Hamiltonian. poj_ 
is orthogonal to the Hamiltonian, and represents the frac- 
tion of the initial spin polarization that does precess. 

III. SPIN-1/2 ELECTRON SYSTEMS 

First we discuss Eq. © for spin-1/2 systems. The 
Hamiltonian describing spin-orbit coupling has the form 
H = (h/2) <j ■ fl{k) and p may be decomposed as p = 
i [n+s(t)-a], where n represents the number density and 
s(t) the spin polarization. Equation ^ has qualitatively 
different solutions depending on the regime under study, 
and they are discussed in detail below. 

A. Exponential decay in the strong momentum 
scattering regime 

A solution to Eq. ^ characterizing relaxation is under- 
stood as exponential decay of the form p(t) = e~ Tst po, 
where T s is generally a second-rank tensor that represents 
the inverse of the spin relaxation time t s . Such a simple 
solution of Eq. ([3]) does not exist in general, but for strong 
momentum scattering (Or p < 1) the RHS of Eq. can 



be neglected. Then substituting for p and H in Eq. ([3]) 
yields the DP expression^ for T s , which may be writ- 
ten as (r s )jj = t p (p, 2 b~ij — f2j57j), where i,j = x,y,z. 
Strong momentum scattering yields exponential spin re- 
laxation and the well-known^ trend r s oc r p _1 . 

B. Oscillations in the ballistic regime 

Previously, most analytical studies have focused on 
strong momentum scattering i 5 ' 6 ' 7 ' 8 ' 13 ' 14 ' 15 ' 16 ' 17 We will 
show that the ballistic and weak momentum scatter- 
ing regimes are far more complex . 37 i 38 ' 39 i 40 In the bal- 
listic limit t p — > oo and Eq. ^ can be solved exactly 
as p(t) = e~ iHt '/ h p 0± e iHt'/n „|_ p which can also be 
obtained from the quantum Liouville equation^ 9 . This 
determines the time evolution of an initial spin polariza- 
tion s(t — 0) = So, i.e., the component of s(t) along 
So. For simplicity So is here assumed independent of k; 
a fc-dependent distribution would not change the results 
qualitatively. From the solution for p in the ballistic limit 
we have 

s(t) ■ s = [1 - (fi • so) 2 ] cos fit + (O ■ s ) 2 , (4) 

where a denotes the unit vector in the direction of a. The 
last term corresponds to p \\ ■ It is best to take a concrete 
example, such as the Hamiltonian of a 2D system on 
a (001) surface with linear Rashba^i and Dresselhaus^ 2 - 
spin-orbit interactions, H = a(o~ x k y — cr y k x ) + (3(cr x k x — 
cTyky). We consider first effective fields fl(k) such that 
the magnitude |Cl(fe)| does not depend on the direc- 
tion of fc, for example either a = or (3 = yields 
|0(fe)| = fl(k). In this case the initial spin polarization 
will simply oscillate with frequency fl. It is helpful to 
visualize a population of spins on the Fermi surface, all 
initially pointing up. If |fi(fe)| is the same at all points fc 
on the Fermi surface, all spins that were in phase initially 
will be in phase again after one precession period. Some 
fraction of the initially oriented spins So, corresponding 
to the last term in Eq. has a nonzero overlap with 
the local field fi(fc) so that the projection of So on fi(fc) 
will be preserved. This fraction is zero if the initial spin 
is out of the plane, but significant if it is in the plane. 

C. Non-exponential decay in the ballistic regime 

The case when |fi(fe)| depends on the direction of k is 
of great relevance to experiment, where spin-orbit cou- 
pling is rarely attributable to a single mechanism. Spins 
on the Fermi surface precess with incommensurable fre- 
quencies and once they are out of phase they never all 
get in phase again (but the polarization fraction due to 
Po|| is conserved.) In our example, if j3 <C ot we can 
write fl Pi fl (1 + 7/ sin 20), where fl = k \J a 2 + (3 2 , the 
small parameter r\ = a(3/[2(a 2 + (3 2 )], and <fi is the po- 
lar angle of fc. The motion of the spins, given by Eq. 
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(HJ, is averaged over the Fermi circle. Consider the term 
cosfH in Eq. Q as an example. The angular average 
yields cos(fii) Jo(rjQ,t), where Jo is a Bessel function of 
the first kind. This function has the form of a decay- 
ing oscillation but it does not reduce to an exponentially 
damped oscillation in any limit. At long times we have 
Jo(r)£lt) — ► y/2/ (TrrjCli) cos(rjClt — 7r/4). A similar, more 
complicated expression in terms of Bessel functions ap- 
plies for the remaining term (Cl ■ sq) 2 cos fit in Eq. ((4]). 

The anisotropy of the Fermi surface introduces a mech- 
anism for non-exponential spin decay^ 3 - with a charac- 
teristic time Trf oc (77O) 1 , referred to as the dephasing 
time Td- For pure spin dephasing, i.e., in the absence of 
momentum scattering, two limiting cases can be distin- 
guished. If fi(fe) _L Sq for all k (e.g., a spin orientation 
perpendicular to the 2D plane on a [001] surface), spin 
dephasing reduces the spin polarization to zero. On the 
other hand, spin dephasing is completely suppressed if 
fi(fc) || So for all k (e.g., a 2D electron system in a sym- 
metric quantum well on a [110] surface with a spin ori- 
entation perpendicular to the 2D planed). In general (in 
particular for 3D systems), an intermediate situation is 
realized where the spin polarization is reduced because 
of dephasing, but it remains finite. The surviving part 
is identified with p || m the initial density matrix. This 
process is referred to as incomplete spin dephasing. 

Analogous results hold for the fc 3 -Dresselhaus modelj 4 ^ 
but the terms leading to dephasing cannot be expressed 
in a simple form due to the complex angular dependence 
of |fi(fe)| on the direction of k. FigureQJa) shows the in- 
complete dephasing of electron spins in bulk GaAs calcu- 
lated using the A: 3 -Dresselhaus model. At long times the 
initial spin polarization settles to a value w 0.33, which 
is independent of any system parameters, including the 
spin-orbit constant. 

D. Weak momentum scattering regime 

In the regime of weak momentum scattering the solu- 
tion to Eq. §3§ may be written approximately as 

p(t) = po|| + e- t/T " e~* m '/ h p al _ e iHt '/ h . (5) 

Since the momentum scattering rate 1/t p is small, the 
term under the overline is taken to lowest order in l/r p . 
The second term on the RHS of Eq. ([5]) describes damped 
oscillations with amplitude decaying exponentially on a 
scale oc t p . This trend is the inverse of that for strong 
momentum scattering and is explained by the following 
argument. If one spin, precessing on the Fermi surface in 
phase with all the other spins, is scattered to a different 
wave vector, it will precess about a different effective field 
and will no longer be in phase with the other spins. Thus 
the combined effect of spin precession and momentum 
scattering — even when the latter is only weak — reduces 
the spin polarization faster. The fraction of the spin po- 
larization corresponding to poll survives. This remaining 



polarization decays via spin-flip scattering (the Elliott- 
Yafet mechanism 9 -^) on much longer time scales. 

An exception occurs when fi(fe) || Sq for all spins. 
This situation is realized, e.g., for a 2D electron system 
in a symmetric quantum well on a [110] surface with a 
spin orientation perpendicular to the 2D plane. For this 
particular case it is well known that spin precession and 
momentum scattering do not affect at all the initial spin 
orientation,^' 22 From the preceding discussion we can un- 
derstand this by noting that the initial density matrix 
po = Po|| commutes with the spin-orbit Hamiltonian and 
Eq. (|3|) shows that p = p || f° r au times. The polarization 
decays eventually via spin-flip scattering.—* 1 ^ 

In the weak momentum scattering regime for non- 
negligible anisotropy the spin decay rate is determined 
by the larger of nCl and l/r p . Momentum scattering in- 
troduces an exponential envelope but in this limit the 
concept of a spin relaxation time is evidently of limited 
use.— In spin- 1/2 systems dephasing will be important 
for high-mobility carriers. Results consistent with our 
findings were obtained experimentally by Brand et 
who studied the oscillatory time evolution of an optically- 
generated spin polarization in a high-mobility 2D elec- 
tron system in a GaAs/AlGaAs quantum well. Similarly, 
in materials in which a nonequilibrium spin density is ex- 
cited, the time evolution of this spin density can be stud- 
ied, for example, by means of magnetic circular dichroism 
techniques^ 

We have assumed an initial spin distribution sharp 
at the Fermi edge. In practice this distribution spans 
a window in fe-space, introducing additional dephasing 
between spins at wave vectors of slightly different magni- 
tudes. In this case even an isotropic spin splitting leads 
to decay, though the polarization due to p Q \\ is still ro- 
bust. For example, in 2D for isotropic spin splitting, 
J dk k cos Qt oc t~ 2 , so that the spin polarization, instead 
of oscillating indefinitely, decays as i~ 2 . 

We have also assumed the initial spin distribution to be 
independent of wave vector. The theory is well-equipped 
to deal with wave-vector dependent spin polarizations. 
(Indeed, the initial spin distribution, contained in the 
density matrix po, is in general wave vector-dependent.) 
The wave vector-independent cases discussed at length 
are intended as examples, and they have been selected as 
more straightforward cases for clarity. 



IV. SPIN-3/2 HOLE SYSTEMS 

Next we discuss spin-3/2 hole systems, which are dif- 
ferent from spin- 1/2 electron systems for several reasons. 
The presence of extra terms in the spin density matrix 
of spin-3/2 systems (in addition to the number density 
and spin polarization) has important consequences for 
spin dynamics * 29 ! 47 Spin-orbit coupling affects the energy 
spectrum in the valence band to a greater extent, and the 
spin orientation often disappears on scales comparable to 
r p . The relation ttr p -C 1 holds less frequently than for 



5 









(a) GaAs: electrons 




— * / 










(b) GaAs: — 


— HH 






- - LH 



1 2 3 5 10 15 20 

time (ns) time (ps) 



FIG. 1: Incomplete spin dephasing of (a) electron spins in 
the fc 3 -Dresselhaus model and (b) of heavy-hole spins (solid 
lines) and light-hole spins (dashed lines) in bulk GaAs in the 
ballistic limit. The vertical axis shows the normalized spin 
polarization s(t) ■ So/|sq|. The initial spin polarization sq 
is assumed to point along [001]. The Fermi wave vector is 
kF = 10 8 m -1 . Note the different time scales in (a) and (b). 



electrons in systems accessible experimentally. 

We consider spin-3/2 holes to be described by the Lut- 
tinger Hamiltonian ) 48 i 49 

H ° = 2^ + ^) ^ " 2 ^ (fc ' 5)2 ] ~ Hc ' (6) 

where mo is the bare electron mass, S the spin oper- 
ator for effective spin 3/2, 7 = (272 + 373 )/5, and 71, 
72, and 73 are Luttinger parameters. He represents the 
anisotropic terms with cubic symmetry 4 ^ which will be 
given below. 

We work first in the spherical approximation in which 
He is neglected. The energy dispersions -Ehh for the 
heavy holes (HHs, spin projection in the direction of k 
is m s — ±3/2) and -Elh for the light holes (LHs, m s — 
±1/2) are £ lh /hhO) = &k 2 (71 ± 2 7 )/(2m ). 



A. Exponential decay in the strong momentum 
scattering regime 

In the strong momentum scattering limit an exponen- 
tial solution pit) — e~ Tst po is possible. The tensor 
r s = t.T 1 1, showing that the relaxation times for all 
spin components are equal, 

1 2 8(h^k 2 \ 2 

where now the frequency fl(k) — (_Elh — Enn)/h = 
2hqk 2 /mo corresponds to the energy difference between 
the HH and LH bands. 2 - Despite the qualitatively dif- 
ferent spin precession, the situation is rather similar to 
electron spin relaxation and can be explained in terms of 
the same random walk picture familiar from the study of 
electron spin relaxation. 



B. Ballistic and weak momentum scattering 
regimes in the spherical approximation 

In the ballistic limit Eq. ([3]) is again solved by p(t) = 
e -iHt>/h p e lHt '/ h . An initial spin polarization will os- 
cillate indefinitely since f2 is the same for all holes on 
the Fermi surface. For weak momentum scattering Eq. 
© applies to holes also. The spin polarization consists 
of damped oscillations, decaying on a time scale oc r p , 
plus a term corresponding to Pqu, which survives at long 
times. po|| does not depend on the Luttinger parameters 
or the Fermi wave vector and will therefore be the same 
in any system described by the Luttinger Hamiltonian. 
This remaining polarization decays via spin-flip scatter- 
ing as discussed in Refs. I34II35I 

C. Cubic-symmetry terms and dephasing 

Dephasing is introduced if the term He with cubic 
symmetry 4 ^ is included in the Luttinger Hamiltonian, 

ft 2 A 

He = {k x ky\^Jx-, Jy \ kyk z {Jy-> J z \ ~r~ k z k x {Jz-i JxY) ? 
mo 

. (8) 

where A = (73 — r y^)/2. The eigenenergies are now given 
by £hh/lh = lih 2 k 2 /2m ± where 

El = 7 fc 4 - 2 7 A [k 4 - 6 (k 2 x k 2 y + k 2 y k 2 + k 2 k 2 x )] + k'A 2 . 

The cubic-symmetry terms contained in Eq. ([5]) are usu- 
ally neglected in charge and spin transport without a 
significant loss of accuracy. However, they play a crucial 
role in spin relaxation in the weak momentum scattering 
regime, which for holes extends over a wide range of k. 

Due to the presence of He, the energy dispersion rela- 
tions and therefore |fi(fe)| depend on the direction of k, 
causing an initial spin polarization to decay even in the 
ballistic limit, where incomplete spin dephasing occurs. 
Our numerical calculations exemplified in Fig.[ljb) show 
that an initial spin polarization falls to a fraction much 
higher than in the electron cases studied. It decays more 
slowly for the LHs, for which the Fermi surface is nearly 
spherical, than for the HHs, for which the Fermi surface 
deviates significantly from a sphere. At long times the 
initial spin polarization settles to a value 0.75, which 
is independent of any system parameters, including the 
Luttinger parameters. 

V. SUMMARY 

In conclusion, we have shown that the decay of spin 
polarization in semiconductors, brought about by the in- 
terplay of spin precession and momentum scattering, de- 
pends strongly on the regime of momentum scattering. 
In the ballistic regime the spin polarization decays via 
a dephasing mechanism which is present due to the fact 
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that the magnitude of the spin-orbit interaction gener- 
ally depends on the direction of the wave vector. This 
mechanism may reduce a spin polarization to zero (com- 
plete dephasing) or a fraction of the initial value (incom- 
plete dcphasing). Weak momentum scattering destroys 
an initial spin polarization, whereas strong momentum 
scattering helps maintain an initial polarization. 
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